void fill0_one(double *u, int n)
//Fills u[1..n]with zeros.
{
  int i;
    for(i=1;i<=n;i++)
      u[i]=0;
}
